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Abstract 

The self-consistent random phase approximation (RPA), based on the framework of relativistic 
energy density functionals, is employed in the study of isovector and isoscalar dipole response in 
^^Ni, ^^^Sn, and ^"^^Pb. The evolution of pygmy dipole states (PDS) in the region of low excitation 
energies is analyzed as a function of the density-dependence of the symmetry energy for a set 
of relativistic effective interactions. The occurrence of PDS is predicted in the response to both 
the isovector and isoscalar dipole operators, and its strength is enhanced with the increase of the 
symmetry energy at saturation and the slope of the symmetry energy. In both channels the PDS 
exhausts a relatively small fraction of the energy-weighted sum rule but a much larger percentage 
of the inverse energy-weighted sum rule. For the isovector dipole operator the reduced transition 
probability B{E1) of the PDS is generally small because of pronounced cancellation of neutron 
and proton partial contributions. The isoscalar reduced transition amplitude is predominantly 
determined by neutron particle-hole configurations, most of which add coherently, and this results 
in a collective response of the PDS to the isoscalar dipole operator. 

PACS numbers: 21.60.Jz, 24.30. Cz, 24.30. Gd, 21.65.Ef, 21.30.Fe 
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I. INTRODUCTION 



Studies of collective excitations play a crucial role in our understanding of complex prop- 
erties of nuclei. The multipole response of unstable nuclei and, in particular, a possible 
occurrence of new exotic modes of excitation in weakly bound nuclear systems presents a 
rapidly growing field of research During the last decade a number of experimen- 

tal studies have been performed on the low-energy electric dipole response in neutron-rich 
medium- heavy and heavy nuclei j^-Q]. The currently available data, however, are not suf- 
ficient to completely determine the nature of observed excitations, and the interpretation of 
the dynamics of low-energy El strength is under discussion. 

The structure and dynamics of low-energy dipole strength, also referred to as pygmy 
dipole state (PDS) or resonance, has extensively been investigated using a variety of the- 
oretical approaches and models Recent studies have made use of the Hartree-Fock 
(HF) plus random phase approximation (RPA) 2lN28|. the self-consistent RPA based on 
the AV18 nucleon-nucleon interaction and phenomenological three-body contact terms 29| . 



the Hartree-Fock-Bogoliubov (HFB) model plus quasiparticle (Q)RPA 
cle RPA plus phonon coupling |37|, the quasiparticle- phonon model 



30 



, quasiparti- 



the quasiparticle 



phonon model including complex configurations of up to three phonons 
second RPA (SRPA) |39|, the relativistic RPA kolW and QRPA 



tic quasiparticle 
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1^ 



the 



44-471, the relativis- 



approach [521. l53|. 



time blocking approximation 48l-l5lj. and the semiclassical Landau- Vlasov 



In addition to the fact that pygmy states present an intrinsically interesting structure 
phenomenon in neutron-rich nuclei, it has been suggested that they might constrain the 
radius of the neutron distribution in medium-heavy and heavy nuclei 9|, and provide in- 
formation about the density dependence of the symmetry energy j4o|. Prompted by recent 
experiments that have reported accurate data on the distribution of PDS in medium heavy- 
nuclei js iSl, and by the first model-independent determination of the neutron-skin thickness 
in ^°^Pb with parity- violating elastic electron scattering 5j|, in the last two years several 
theoretical studies have explored the relation between the PDS, the density dependence of 
the symmetry energy, and the neutron skin 23|, |2J, |28|, |43|] . These studies are based on the 



microscopic framework of nuclear energy density functionals (EDF) plus the (Q)RPA. 
Starting from a representative set of Skyrme effective forces and meson-exchange effective 
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Lagrangians, Carbone et al. 23|] performed an RPA analysis of the correlation between the 
density dependence of the nuclear symmetry energy, the neutron skin, and the percentage 
of energy-weighted sum rule (EWSR) exhausted by the PDS in ^^Ni and ^^^Sn. In com- 
parison with available data it was possible to constrain the value of the derivative of the 
symmetry energy at saturation, and use this constraint to determine the neutron-skin radii 
of ^^Ni, ^^^Sn, and ^°^Pb. In contrast to this result, Reinhard and Nazarewicz 2^, using the 



covariance analysis to identify observables and pseudo-observables that correlate with the 
neutron skin, suggested that the neutron skin of ^°'^Pb is strongly correlated with the dipole 
polarizability but very weakly correlated with the low-energy electric dipole strength. This 
finding has recently been challenged by Piekarewicz 43| in an analysis of the distribution of 
electric dipole strength in ^^Ni using a relativistic RPA with a set of effective interactions 
that predict significantly different values for the neutron-skin thickness of ^°^Pb. The results 
suggest a strong correlation between the dipole polarizability of ^^Ni and the neutron-skin 
thickness of ^°^Pb, but also a correlation just as strong between the neutron-skin thickness of 
^°^Pb and the fraction of the dipole polarizability exhausted by the pygmy dipole strength. 
In a very recent study performed using the self-consistent Skyrme HP plus RPA approach, 
X. Roca-Maza et al. 28| have analyzed the isospin character, the degree of collectivity, 
and the sensitivity to the slope of the nuclear symmetry energy, of the low-energy dipole 
response in ^^Ni, ^^^Sn, and ^'^^Pb. It has been shown that both the isoscalar and isovector 
strength functions display a low-energy peak that is enhanced and shifted to higher excita- 
tion energies with increasing values of the slope of the symmetry energy at saturation. The 
degree of collectivity associated with the RPA state(s) contributing to this peak differs in 
the isoscalar and isovector channels. Much more collectivity in the PDS is predicted in the 
response to the isoscalar dipole operator. 



In this work we perform an analysis similar to that of Ref. [28[, but using a more sys- 
tematic set of effective nuclear interactions. Namely, to analyze the model dependence of 
the predicted PDS, X. Roca-Maza et al. |28| employed three different Skyrme parameter 
sets: SGII, SLy5 and SkI3. These interactions span a broad range of values of the slope 
of the nuclear symmetry energy at saturation, but they also differ in other characteristics 
in a nonsystematic way. A consistent set of effective interactions was used by Piekarewicz 



in Ref. 
channel. 



43| to analyze the distribution of PDS, but only for ^^Ni and only in the isovector 
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II. THEORETICAL FRAMEWORK 



The present analysis employs the fully self-consistent relativistic random phase approx- 
imation (RRPA) based on the framework of relativistic energy density functionals [55]. In 
the relativistic mean-field (RMF) + RPA model effective interactions are implemented in a 
fully consistent way: effective Lagrangians with density- dependent meson-nucleon couplings 



are employed 



56 



57| and the same interactions are used both in the RMF equations that 



determine the ground state, and in the matrix equations of the RRPA. The full set of RRPA 
equations is solved by diagonalization. The results are excitation energies Ex and the corre- 
sponding forward- and backward-going amplitudes, and Y^, respectively, that are used 
to evaluate the reduced transition probabihty from an excited state \JX) to the ground state: 

1 



(EJ) 



2Ji + l 



(1) 



where /i and /x' denote single-nucleon states. Discrete spectra are averaged with a Lorentzian 
distribution of arbitrary width (1.5 MeV in the present calculation). The electric El response 
is calculated for the isovector dipole operator: 
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and the isoscalar dipole operator: 
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The inclusion of the second term in the isoscalar operator Eq. ([3]), with rj = 5(r^)/3, en- 
sures that the corresponding strength distribution does not contain spurious components 



associated to the center-of-mass motion 



58| . The strength function reads 



S{E) = Y.\MQ^jm\'S{E-E^). 



(4) 



El, is the energy of the RPA state |z/), and the moments of the strength distribution are 
defined as 

[dEE'S{E)=j:Etmm\0)\'. (5) 



rrik 
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In the following we analyze the occurrence and structure of the PDS in the isovector and 
isoscalar dipole response of ^^Ni, ^^^Sn, and ^"^Pb, in relation with the density dependence 
of the symmetry energy. In linear order with respect to the nuclear matter density p, the 
symmetry energy S{p) is determined by its value at saturation density S{po) = a^, and by 
the derivative at saturation density: 



5'(p)U,^ — , (6) 
and this relation defines the "slope" parameter L. Using data on the percentage of the 



EWSR associated with the PDS in ^^Ni 



15| and ^^^Sn {q], Carbone et al. 23| constrained 



the value of the slope parameter: L = 64.8 ±15.7 MeV, in accordance with values previously 
determined with different types of analyses and/or other methods based on nuclear struc- 
ture and heavy-ion experiments. Here we employ the framework of relativistic energy density 
functionals represented by effective Lagrangians with density-dependent meson-nucleon ver- 
tex functions. The parameters of the very successful effective interactions DD-MEl 56[ and 



DD-ME2 



57| . in particular, were adjusted simultaneously to empirical properties of sym- 



metric and asymmetric nuclear matter, and to binding energies and charge radii of twelve 
spherical nuclei. Data on excitation energies of isoscalar monopole (ISGMR) and isovector 
dipole giant resonances (IVGDR) were also used to determine the compressibility modulus 
and asymmetry energy at saturation, as well as available data on differences between neutron 
and proton radii in Sn isotopes and ^''^Pb. Numerous calculations have shown that these 
interactions provide accurate results for ground-state properties of spherical and deformed 
nuclei, as well as excitation energies of giant resonances. Pertinent to the present analysis, 
the relativistic RPA with the DD-ME2 effective interaction predicts the dipole polarizability: 

Stt o 



aD 



9 



■e m_i 



(7) 



(directly proportional to the inverse energy-weighted moment for ^°^Pb: 20.8 fm^, in 

very good agreement with the recently measured value: ao = (20.1 ± 0.6) fm^ 20 |. 



III. RESULTS 

In an earlier study [59*] we used the relativistic RPA with density-dependent meson- 
nucleon effective interactions to provide a microscopic estimate of the nuclear matter com- 



5 



> 

0) 



E 
> 

CQ 




10 15 20 25 30 35 40 

E (MeV) 



20 30 

E (MeV) 



40 50 



FIG. 1: (Color online) The relativistic RPA isovector (left column) and isocalar (right column) 
dipole strength distributions in ®^Ni, ^'^^Sn, and ^"^Pb, calculated with a set of five effective inter- 
actions that differ in the value of the symmetry energy at saturation (04) and the corresponding 
slope parameter L. 

pressibility and symmetry energy in relativistic mean-field models. Starting from the param- 



56| , three families of interactions with nuclear matter incompressibility 



eterization DD-MEl 

i^nm = 230, 250, and 270 MeV were constructed. For each value of five interactions 
were adjusted with = 30, 32, 34, 36 and 38 MeV, respectively. For particular values of 
-f^nm and 04, the remaining parameters of these interactions were fine tuned to empirical 
properties of nuclear matter, and to the binding energies and charge radii of ten spherical 
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nuclei. It was noted that, to reproduce the empirical masses oi N Z nuclei, larger values 
of 04 necessitate an increase of the slope parameter L of the symmetry energy. These effec- 
tive interactions were used in Ref. [9|] to study a possible correlation between the observed 
PDS in ^^"'^'^^Sn and the corresponding values for the neutron skin thickness. In addition to 
DD-ME2 (i^nm = 251 MeV), the set of effective interactions with K^^ = 250 MeV and 04 = 
30, 32, 34, 36 and 38 MeV will here be used to study PDS in the isovector and isoscalar 
dipole response of spherical closed-shell nuclei. 

For this set of effective interactions, in Fig. [1] we display the total isovector (IV) and 
isoscalar (IS) RRPA dipole strength functions for ^'^Ni, ^^^Sn, and ^°^Pb. The IV dipole 
response in all three nuclei is, of course, dominated by the collective giant dipole resonance 

GDR) peaked in the high-energy region around 15 MeV. As it was already shown in Refs. 

55|,|60|, as a result of the increase of L with 04 the excitation energy of the IV GDR decreases 
with increasing S'(po) = cla, because this increase implies a reduction of the symmetry energy 
at low densities characteristic for surface modes. In addition, one notices an enhancement 
of El strength in the low-energy region below 10 MeV. This pygmy dipole strength (PDS) is 
very sensitive to the density dependence of the symmetry energy, and is strongly enhanced 



by the increase of L and 04. As shown in Ref. [28|], the same pygmy states are also present 
in the IS strength functions, and this is a clear evidence of the mixed isovector-isoscalar 
nature of the PDS. The isoscalar El strength distributions display a characteristic bimodal 
structure with two broad components: one in the low-energy region close to the IV GDR 
(^ 2huj), and the other at higher energy close to the electric octupole resonance Shu). 
Theoretical analyses have shown that the high-energy component represents compressional 
vibrations 



58 



61 



62| . The high-energy IS GDR is a second order effect, built on 3hu! 



or higher configurations, and corresponds to a compression wave traveling back and forth 
through the nucleus along a definite direction. Some states comprising the broad structure 
in the low-ener gy r egion correspond to vortical nuclear flow associated with the toroidal 



dipole moment |52|, I63l-l65|. However, as pointed out in a study of the interplay between 
compressional and vortical nuclear currents 6^, a strong mixing between compressional and 
vorticity vibrations in the isoscalar El states can be expected up to the highest excitation 
energies in the region ^ Shu. Finally, the lowest peaks in the isoscalar strength functions 
are associated to the PDS. A very interesting result, also pointed out in Ref. is that 
relatively to the corresponding GDR, the PDS is much more pronounced in the isoscalar 
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channel. In fact, for ^^Ni the pygmy state has the largest B{E1) value among the isoscalar 
states, and the IS strength function, except the structure around 15 MeV, is not sensitive to 
the variation of the density dependence of the symmetry energy. One expects more isospin 
mixing in the two heavier nuclei (see also the analysis of Ref. 45|]), and this is reflected in 
the enhancement of the B{E1) values with the increase of L and 04, for all states in the 
response of ^^^Sn and ^°^Pb to the isoscalar dipole operator Eq. ([3]). This enhancement is 
more pronounced below the IS GDR, that is, in the region around 15 MeV where one expects 
mixing with the IV GDR, and for the PDS. 

The isoscalar dipole strength functions for ^°^Pb, shown in Fig. |21 have been calculated 
for a complementary set of interactions: with the isovector parameters fixed at 04 = 32 
MeV and L = 46.5 MeV, the four interactions differ only in the value of the nuclear matter 
incompressibility: K^^a = 210, 230, 250, and 270 MeV. The high-energy component is 
obviously very sensitive to the choice of the nuclear matter incompressibility. Fragmentation 
in the region 20 - 30 MeV is somewhat reduced with increasing -ft"nm, and the centroid of this 
structure is shifted to higher energies. Varying -R'nm has very little effect on the pygmy dipole 
state: the energy does not change, and only the B{E1) is slightly reduced with increasing 
incompressibility as the total strength gets shifted to higher energy. From this result one 
deduces that the admixture of 3hu compressional mode in the isoscalar PDS is indeed very 
small. An intermediate situation is found for the pronounced resonant structure in the energy 
region between 10 and 15 MeV. As a result of the stronger mixing with the compressional 
mode the centroid is shifted to higher energy, but at the same time the transition strength 
is reduced considerably with increasing K^m- 

Even though the set of effective interactions that we use here spans a relatively broad 
interval of values of the parameters 04 and L that characterize the density dependence of the 
symmetry energy, the effect on the excitation energy of the IV GDR is not large, especially 
in a heavy nucleus like ^°^Pb. In fact, as shown in Fig. [H the IV GDR peak energy in ^osp]^ jg 
lowered by less than 1 MeV in the interval between 04 = 30 MeV, L = 30 MeV and 04 = 38 
MeV, L = 111 MeV. A better distinction between these effective interactions is obtained by 
considering predicted values for the electric dipole polarizability of ^osp]^ ([7])), in 

comparison with the experimental value an = (20.1 ±0.6) fm^ 20|]. The results are shown in 
Fig. [3l In addition to the values obtained with the five effective interactions 04 = 30, 32, 34, 
36, and 38 MeV, we also include the polarizability calculated with the interaction DD-ME2 
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FIG. 2: (Color online) The relativistic RPA isoscalar dipole transition strength in ^'^^Pb. The 
four curves correspond to effective interactions that differ only in the value of the nuclear matter 
incompressibility. 



(04 = 32.3 MeV and L = 51.5 MeV) that was adjusted independently in Ref. 57|. The 
result closest to the experimental value corresponds to the effective interaction with 04 = 32 
MeV and L = 46.5 MeV, and the value predicted by DD-ME2: an = 20.8 MeV is just 100 
keV outside the experimental error bar. The slope parameters of these two interactions are 

I deduced from the 



slightly below the value L = 64.8 ± 15.7 MeV that Carbone et al. 

percentage of the EWSR associated with the PDS in ^^Ni |l^ and ^^^Sn [9|]. From Fig. [3] 
it is apparent that only the values of predicted by the two interactions with 04 = 36 
MeV, L = 85.5 MeV and 04 = 38 MeV, L = 110.8 MeV are in serious disagreement with 
experiment. 

For ^^Ni in Fig. H] we display the moments of the dipole strength distributions mo, mi and 
m_i for the isovector (panel (a)-(f)) and isoscalar (panel (g)-(i)) channels. E = 12 MeV is 
the mo del- dependent, although obvious choice for the excitation energy that separates the 
low-energy PDS region from the energy interval dominated by the collective giant resonance 
(cf. Fig.[l]). In the left column we plot the absolute values of the moments of the isovector and 
isoscalar strength distributions in the PDS region, whereas the column on the right displays 
the percentages of the total moment exhausted in the PDS interval, as a function of the 
slope parameter L of the six effective interactions used in the present analysis. The shaded 
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FIG. 3: Predictions for the electric dipole polarizability Eq. ([7]) of ^''^Pb, calculated with a set of 
five effective interactions that differ in the value of the symmetry energy at saturation (04) and the 



57|. The 



corresponding slope parameter L (circles), and the effective interaction DD-ME2 (star) 
shaded area delineates the experimental constraint = (20.1 it 0.6) fm'^ [20]. 



area denotes the experimental result and the corresponding; uncertainty of the percentage 
of EWSR associated with the PDS in ^SNi 
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23| . We note that the absolute values of 



moments in the low-energy region, as well as the percentages of the total moment exhausted 
in the PDS interval, increase linearly with the slope of the symmetry energy. In both the 
isovector and isoscalar channels the smallest percentage exhausted by the PDS is for the 
mi moment. This is because mi is energy-weighted and, therefore, dominated by the main 
giant-resonance structure in the high-energy region. Much more sensitive to the PDS is the 
moment m_i, directly proportional to the total polarizability Eq. ([7]), because the inverse 
energy- weighting enhances the low-energy part of the response. In the isovector channel, 
in particular, the percentage of the total polarizability exhausted by the PDS is more than 
three times the percentage of the EWSR exhausted by the low-energy pygmy strength. 
These results are in agreement with those reported by Piekarewicz [43,], obtained using a 
different class of relativistic effective interactions but with a similar range of values for the 
slope parameter of the symmetry energy. What is new here are the results for the moments 
of the isoscalar strength distributions shown in the lower panel of Fig. HI In this case the 
PDS exhausts more than 10% of the EWSR, and more than 50% of the m_i moment, for all 
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FIG. 4: Moments mo, mi and m_i of the isovector (panel (a)-(f)) and isoscalar (panel (g)-(i)) 
dipole strength distributions of ^^Ni. Absolute values of the moments in the low-energy PDS 
region below E < 12 MeV (left column), and the percentages of the total moments exhausted by 
the PDS (right column), as a function of the slope parameter L of the six effective interactions 
used in this work . 
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values of the slope parameter L, in agreement with the results of X. Roca-Maza et al. [28 1 
obtained using Skyrme effective interactions. Similar results are also found for the moments 
of the isovector and isoscalar dipole strength distributions of ^'^^Sn and ^°'^Pb. 

An important issue is the degree of collectivity of the PDS. In most of the early theoretical 
studies the measure of collectivity was associated with the number of particle-hole config- 
urations that significantly contribute to the RPA amplitude of the principal pygmy state. 

Such a measure, however, does not take into account the coherence of these contributions. 

I I] 

Lanza et al. [22j] analyzed the collectivity of PDS by considering the partial contributions to 
the reduced transition amplitude. Namely, the reduced transition probability Eq. ([T]) from 
the ground state to the excited state | JA) can be written as: 



B'\EJ) 



ph 



(8) 



where the summation is over all particle-hole configurations that build the excited RPA state. 
This expression shows that the reduced transition probability is determined by the number 
of configurations that contribute with a significant weight, as well as by the coherency 
(relative sign) of these contributions. The RPA state can be considered as a resonance 
if the corresponding reduced transition amplitude is composed of more than just a few 
particle-hole partial amplitudes similar in magnitude and with the same relative sign. By 
comparing the contributions of the partial amplitudes A^r^'^ jEl) for the pygmy state and 



22I ] concluded that, although the 



the isovector giant dipole resonance in ^^^Sn, Lanza et al. 
main pygmy state can indeed be formed by many particle-hole configurations, the reduced 
transition probability B{E1) is generally small because the corresponding partial amplitudes 
cancel out to a large extent. This is in contrast to the GDR for which many configurations 
add coherently to build a large transition amplitude. In their study of low-lying dipole 
response in ^^Ni, ^^^Sn, and ^°'^Pb, X. Roca-Maza et al. [2^ have found that the isovector 
reduced amplitude of the PDS is characterized by destructive interference of a relatively 
small number of different particle-hole configurations, and the resulting reduced transition 
probability does not exceed ^2 — 4 single-particle units. In the isoscalar channel the largest 
contributions to the reduced transition amplitude predominantly arise from neutron particle- 
hole configurations, and most of them add coherently. This results in a collective response of 
the PDS to the isoscalar dipole operator, characterized by a reduced transition probability 
of ~ 10 — 20 single-particle units, for all three nuclei and the three interactions that were 
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used in the analysis of Ref. [28 1. 

The difference between the PDS in the isovector and isoscalar dipole response of ^^^Sn is 
illustrated in Fig. O where we plot the largest partial contributions of neutron and proton 
particle-hole configurations to the isovector (left) and isoscalar (right) reduced transition 
amplitude of the main pygmy state at 7.81 MeV excitation energy, calculated using the 
effective interaction with 04 = 32 MeV and L = 46.5 MeV. The partial reduced transition 
amplitudes in units of fm (isovector) and fm^ (isoscalar), are plotted as a function of the 
unperturbed energy of the corresponding particle-hole configurations. Only amplitudes with 
a magnitude larger than ^ 0.01 fm (isovector), and ^0.1 fm^ (isoscalar) are shown in the 
figure. About 7 to 8 neutron particle-hole configurations display a non-negligible partial 
transition amplitude, whereas only two proton configurations contribute significantly to the 
isovector and isoscalar amplitudes. One notices, however, the destructive interference of pro- 
ton and neutron particle-hole configurations in the isovector channel (cf. also Table |T]). The 
dominant neutron amplitudes in the PDS energy region correspond to the configurations 
3si/2 3p3/2, 2^3/2 3pi/2, 3si/2 3^1/2, and 2^3/2 3p3/2, all with positive sign. The 
two large negative neutron amplitudes in the GDR region that correspond to the configura- 
tions l/iii/2 — 1^13/2 and 1(77/2 1^9/2, are almost exactly canceled by the positive proton 
contributions from 1(79/2 l/iii/2 and I/5/2 ^97/2, respectively. It is interesting to note 
that, as shown in Fig. [6l it is precisely the partial transition amplitudes of the latter two 
neutron and two proton configurations that add coherently (here with a negative sign) to 
produce the large collective B(E1) of the IVGDR state at 15.24 MeV. The repulsive residual 
interaction, of course, gathers these contributions and shifts them to higher energy to build 
the collective GDR. Therefore, while the large proton and neutron reduced amplitudes in 
the high-energy region add coherently to build the GDR, for the PDS state at 7.81 MeV the 
same particle-hole configurations with positive proton and negative neutron contributions 
cancel each other, so that the rather small reduced transition amplitude of the PDS state 
is dominated by just a few neutron particle- hole configurations (note that the unperturbed 
PDS strength is shifted to lower energy by the residual interaction). 

The picture is very different in the isoscalar channel (cf. also Table HI])- Also in this 
case the dominant neutron configurations for the PDS are 3si/2 -> 3p3/2, 2(^3/2 — )■ 3pi/2 
3si/2 — )■ 3pi/2, and 2^3/2 — ?■ 3p3/2 (positive sign for the partial transition amplitude), but 
here the contributions from 1/111/2 1^13/2 and 1(^7/2 1/^9/2 are also positive even though 
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FIG. 5: (Color online) Partial contributions of neutron and proton configurations to the isovector 
(left) and isoscalar (right) reduced transition amplitude of the PDS state at 7.81 MeV excitation 
energy in ^^^Sn for the effective interaction with = 32 MeV and Xnm = 250 MeV, as functions 
of the unperturbed energy of the particle-hole configurations. Only amplitudes with a magnitude 
larger than « 0.01 fm (isovector), and ~ 0.1 fm^ (isoscalar) are shown in the figure. 

an order of magnitude smaller. Because most neutron configurations add coherently to the 
reduced transition amplitude, the total neutron amplitude ~ 84 fm^ is much larger than 
the corresponding proton contribution ^17 fm^. In the response to the isoscalar dipole 
operator all major neutron and proton partial transition amplitudes for the PDS state at 
7.81 MeV are of the same sign, and this coherence leads to the large B{E1) value shown in 

Fig.m 

We have verified that similar results are also obtained with the other effective interactions 
used in this work. Tables [T] and [III hst the dominant partial reduced transition amplitudes 
for the PDS in the isovector and isoscalar channel in i32Sn, respectively, calculated for the 
five effective interactions with = 30, 32, 34, 36, and 38 MeV and the corresponding 
increasing values of the slope parameter L. In the isovector channel both the neutron and 
proton reduced transition amplitudes increase with 04 and L (cf. Fig. [1]). Table [T] shows that 
the contributions of the neutron and proton particle-hole configurations increase by 72% 
and 70%, respectively from 04 = 30 MeV, L = 30 MeV to 04 = 38 MeV, L = 110.8 MeV. 
The increase of 04 and L simply corresponds to an enhancement of the restoring force for 
isovector oscillations. However, because of the partial cancellation of proton and neutron 
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FIG. 6: (Color online) Partial contributions of neutron and proton configurations to the isovector 
GDR state at 15.24 MeV excitation energy in ^'^^Sn for the effective interaction with 04 = 32 MeV 
and Knin = 250 MeV, as functions of the unperturbed energy of the particle-hole configurations. 
Only amplitudes with a magnitude larger than ~ 0.01 fm are shown in the figure. 

transition amplitudes, the isovector B{E1) of the PDS is generally small for all values of 
04 and L. The occurrence and evolution of strength of the PDS has been associated with 
the neutron excess. The correlation between neutron-skin thickness and the parameters 04 
and L has been studied using a variety of non-relativistic and relativistic mean-field models 







66|, 



671]. It has been shown that for a given nucleus the thickness of the neutron skin 



increases linearly with 04 and L. In fact, using the same set of effective interactions, in 



Ref. [9] we have calculated the increase of the neutron skin in ^'^^Sn from r„ 



0.19 fm 



for 04 = 30 MeV, L = 30 MeV to r„ - rp = 0.36 fm for = 38 MeV, L = 110.8 MeV. A 
number of studies, including the present, have shown that the strength of the pygmy state 
in the isovector response is enhanced with increasing 04 and L, and this enhancement has 
been correlated with the corresponding increase of the neutron-skin thickness. 

The situation is different in the isoscalar response. In this case the increase of the total 
proton transition amplitude between 04 = 30 MeV, L = 30 MeV and = 38 MeV, L = 
110.8 MeV (39%) is much larger than that of the neutron transition amplitude (18%). 
However, the total contribution of proton configurations to the overall amplitude is small, 
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TABLE I: Partial neutron and proton isovector reduced transition amplitudes A^j^ ^{El) for the 
main pygmy state in I32sn. Included in the table are only amplitudes with a magnitude larger 
than about 0.1 fm, calculated for the five effective interactions with 04 = 30, 32, 34, 36, and 38 
MeV, and the correspondingly increasing values of the slope parameter L. The total neutron (n) 
and proton (p) amplitudes are obtained by summing over all configurations, including also those 
not listed in the table. 



04 (MeV) 






30 


32 


34 


36 


38 


E (MeV) 






7.94 


7.81 


7.71 


7.67 


7.62 




config. 








neutron 


3Si/2 


3p3/2 


0.72 


0.72 


0.70 


0.70 


0.69 




l%l/2 - 


' l«13/2 


-0.59 


-0.60 


-0.60 


-0.62 


-0.65 




24/2 ^ 


3pi/2 


0.20 


0.22 


0.23 


0.24 


0.25 




I97/2 


1/19/2 


-0.19 


-0.20 


-0.19 


-0.20 


-0.20 




3si/2 


3pi/2 


0.18 


0.17 


0.16 


0.16 


0.16 




24/2 ^ 


3p3/2 


0.11 


0.12 


0.13 


0.14 


0.14 




24/2 ^ 


2/7/2 


-0.12 


-0.11 


-0.10 


-0.08 


-0.05 


n-total 






0.32 


0.38 


0.42 


0.49 


0.55 


proton 


159/2 


1^11/2 


0.49 


0.53 


0.55 


0.60 


0.66 




1/5/2 ^ 


157/2 


0.18 


0.20 


0.20 


0.23 


0.25 


p-total 






0.53 


0.64 


0.70 


0.81 


0.91 



varying from 16% for 04 = 30 MeV, L = 30 MeV, to 19% for 04 = 38 MeV, L = 110.8 
MeV. The structure of the PDS in the isoscalar response is, therefore, dominated by a 
relatively large number of neutron particle- hole configurations that, together with proton 
configurations, add coherently to build a reduced transition probability of collective strength. 
We also notice that with the systematic increase of 04 and L the excitation energy of the 
PDS actually decreases, in contrast to the result reported in Ref. [28| but consistent with 
the behavior of the corresponding IV GDR (cf. Fig. [1]). 

The discrete isovector and isoscalar B{E1) spectra for ^'^^Sn are compared in Fig. [71 
for the effective interaction with 04 = 32 MeV, L = 46.5 MeV, and ifnm = 250 MeV. 
At energies below ^ 20 MeV many RRPA states display pronounced B{E1) values both 
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TABLE II: Same as described in the caption to Table [H but for the partial isoscalar reduced 
transition amplitudes AEj^^{El) with a magnitude larger than about 1 fm^. 



OA (MeV) 

Uj/^ \ i-V J-^i; V ( 






30 


32 


34 


36 38 


E (MeV) 






7.94 


7.81 


7.71 


7.67 7.62 




config. 






TlPn tTOTl 


3si/2 


3P3/2 


33.43 


32.96 32.45 31.75 31.18 




2^3/2 ^ 


3Pl/2 


12.85 


14.13 15.08 15.98 16.48 




3si/2 


3pi/2 


Q 1 n 


8.72 


8.11 


7.87 7.75 




2^3/2 -> 


3p3/2 


7.11 


7.73 


8.41 


8.75 8.78 




157/2 


lhg/2 


3.86 


3.72 


3.51 


3.46 3.43 




1^11/2 - 


' 1^3/2 


3.24 


2.88 


2.45 


2.27 2.18 




2^5/2 ^ 


3p3/2 


1.92 


2.60 


3.07 


3.78 4.34 




1^11/2 - 


^ 3ii3/2 


lis 

i.iO 


1.32 


1.41 


1.57 1.69 




2^5/2 ^ 


2/7/2 


-1.16 


-1.10 


-1.11 


-0.86 -0.57 




1^7/2 ^ 


2/5/2 


0.84 


1.12 


1.33 


1.77 2.23 




1^11/2 - 


> 2^9/2 


0.61 


1.08 


1.45 


2.06 2.60 




1^11/2 - 


^ 2ii3/2 


0.81 


0.91 


0.98 


1.08 1.15 


n-total 






79.70 


83.89 86.29 90.76 94.31 


proton 


I59/2 


1^11/2 


5.36 


6.19 


6.68 


7.73 8.74 




1/5/2 -> 


I57/2 


3.32 


3.82 


4.11 


4.76 5.38 


p-total 






15.69 16.95 17.83 19.89 21.77 



in the response to the isovector and isocalar operators. This clearly points to a strong 
mixing between the isovector and isoscalar channels, not only for the PDS but also at higher 
energies. Several states obviously manifest predominantly isoscalar or isovector character, 
e.g. the pygmy state at 7.81 MeV is mainly isoscalar, whereas the main GDR state at 15.24 
MeV is primarily isovector. At excitation energies above ~ 20 MeV the B{E1) spectra 
basically contain only isoscalar transitions that form the broad resonance corresponding to 
the dipole compressional mode. The present results are in agreement with those obtained in 
a previous RQRPA study of PDS based on transition densities 45|]. Theoretical studies have 
also confirmed experimental evidence for the splitting of low-energy El strength into two 
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FIG. 7: (Color online) Discrete B{E1) spectra of the isovector and isoscalar dipole operators in 
^^^Sn, calculated using the effective interaction with 04 = 32 MeV (L = 46.5 MeV), and i^nm = 250 
MeV. 

groups of of states with different isospin structure, one that is excited in (a, a'7) and (7, 7O 
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reactions, and another group of states that is excited only in (7, 7') reactions 0, [l 
However, as the present analysis indicates, rather than finding only two groups of states in 
the low-energy region, one could actually expect a more complex pattern of mixed isoscalar 
and isovector states across the dipole excitation spectra up to ~ 20 MeV. 

Finally, to illustrate the role of proton-to-neutron asymmetry in the evolution of PDS in 
the isoscalar channel, in Fig. |8] we display the isoscalar dipole response for the even-even 



isotopes 



50-68 



Ni. The calculation has been performed with the RQRPA based on the rel- 



ativistic Hartree-Bogoliubov (RHB) model, with the DD-ME2 effective interaction in the 
particle-hole channel, and the Gogny DIS force in the particle-particle channel. For each 
nucleus the RQRPA strength distribution is shown in comparison with the unperturbed 
RHB response. We find that for N ^ Z nuclei the transition strength is strongly frag- 
mented, whereas in other isotopes, either with proton or neutron excess, more pronounced 
peaks are calculated both in the low-energy region and at higher energies. The low-energy 
strength is very sensitive to the proton-to-neutron asymmetry. The reason is its underlying 
nature, i.e. the transition densities for ^'^■^^Ni are dominated by proton configurations, and 
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FIG. 8: (Color online) Evolution of the isoscalar dipole strength distribution in the even-even 



isotopes 



50-68 



Ni, calculated with the RHB+RQRPA model. The DD-ME2 effective interaction is 



used in the particle-hole channel, and the Gogny DIS force in the particle-particle channel. The 
dot-dashed (red) curve corresponds the unperturbed RHB response. 

their contribution decreases towards N ^ Z. On the other hand, for ^^~^^Ni the PDS is 
predominantly built from neutron configurations. Also the high-energy strength around 30 
MeV that corresponds to the Shu compressional mode appears to be rather sensitive to 
the neutron excess. In isotopes beyond ^°Ni a pronounced peak develops at an excitation 
energy of ~ 31 MeV. In all isotopes the structure between the PDS and the high-energy 
compressional peak essentially coincides with the unperturbed RHB response. 



IV. SUMMARY 



In conclusion, we have employed the fully self-consistent random phase approximation 
(RPA), based on the framework of relativistic energy density functionals, to analyze the 
isovector and isoscalar dipole response in ^^Ni, ^^^Sn, and ^''^Pb. In particular, the evolution 
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of pygmy dipole states (PDS) in the region of low excitation energies has been analyzed 
as a function of the density-dependence of the symmetry energy for a set of relativistic 
effective interactions. These interactions, adjusted to empirical properties of nuclear matter, 
binding energies and charge radii of ten spherical nuclei, principally differ in their isovector 
properties. They span a broad range of values of the two parameters that determine the 
density- dependence of the symmetry energy in nuclear matter: the symmetry energy at 
saturation density a^, and the "slope" parameter L proportional to the first derivative of the 
symmetry energy at saturation. The present study has confirmed results recently obtained in 
the framework of non-relativistic and relativistic mean-field plus RPA models: (i) the range 
of values of the slope parameter L constrained by the measured electric dipole polarizability 
of 208p]2 |20:j is consistent with the values deduced from the percentage of EWSR associated 
with the PDS in ^^Ni [isl and ^^^Sn |9]; (ii) the occurrence of PDS is predicted in the 
response to both the isovector and isoscalar dipole operators, and its strength is enhanced 
with the increase of 04 and L; (iii) in both channels the PDS exhausts a relatively small 
fraction of the EWSR but a much larger percentage of the inverse energy-weighted sum rule 
(m_i moment of the strength distribution): ^ 20% in the isovector channel and more than 
50% in the isoscalar channel; (iv) for the isovector dipole operator the reduced transition 
probability B{E1) of the PDS is generally small because of pronounced cancellation of 
proton and neutron particle-hole contributions to the reduced transition amplitude; (v) the 
isoscalar reduced transition amplitude is predominantly determined by neutron particle-hole 
configurations, most of them add coherently and, together with proton contributions, this 
results in a collective response of the PDS to the isoscalar dipole operator. The study of 
performed with the Skyrme HF+RPA model, as well as the present analysis based 



Ref. 



28 



on relativistic effective interactions, have shown that the PDS is much more pronounced 
in the response to the isoscalar dipole operator, and can be considered as a resonance 
state only in the isoscalar channel. This result indicates that isoscalar probes might be a 
more appropriate tool for experimental studies of pygmy dipole states in neutron-rich nuclei 
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